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ABSTRACT 

We report a further analysis of the ratio of deuterium to hydrogen (D/H) using HST 
spectra of the z = 0.701 Lyman limit system towards the QSO PG1718+481. Initial 
analyses of this absorber found it gave a high D/H value, 1.8-3.1 x 10""^ ( Webb et aj] 
Il998j) . inconsistent with several higher redshift measurements. It is thus important 
to critically examine this measurement. By analysing the velocity widths of the Di, 
Hi and metal lines present in this system, Kirkman ct al. (2001) report that the 
additional absorption in the blue wing of the Lya line can not be D i, with a confidence 
level of 98%. Here we present a more detailed analysis, taking into account possible 
wavelength shifts between the three sets of HST spectra used in the a nalysis. We find 
that th e constraints on this system are not as strong as those claimed bv lKirkman et all 
l)2001() . The discrepancy between the parameters of the blue wing absorption and the 
par ameters expe c ted fo r D i is marginally worse than la. 

iTvtler et all (|1999D commented on the first analysis of lWebb et all ()1997l IT99^ . 
reporting the presence of a contaminating lower redshift Lyman limit system, with 
log[A'"(Hi)] = 16.7 at z = 0.602, which biases the iV(Hi) estimate for the main system. 
Here we show that this absorber actually has log[Af(Hi)] < 14.6 and does not impact 
on the estimate of Af(Hi) in the system of interest at z = 0.701. 

The purpose of the present paper is to highlight important aspects of the analysis 
which were not explored in previous studies, and hence help refine the methods used 
in future analyses of D /H in quasar spectra. 

Key viTords: 



1 INTRODUCTION 

Virtually all the deuterium (D) and the vast majority of H 
and ''He we observe today were produced during big bang 
nucleosynthesis (BBN). BBN also produced small amounts 
of ^He and ^Li. In standard BBN theory, the primordial 
abundances of all these light elements depends solely on 
the cosmological baryon density, Qi,. This allows Qijh^ to 
be measured directly by finding the primordial ratio of any 
two of these light elements^ . The ratio D /H is very sensitive 
to Qbh'^ , and of all the BBN light elements, D/H potentially 
gives the best constraints on fi^. 

QSO absorption systems provide a unique way 
to mea s ure the primordial abundance of deuterium 
jAdamj Il97d) . Theoretically there are processes other 
than BBN that create D, but it has been shown 
iEpstein. Lattimer fc Schramml fl974l) that the D produc- 
tion due to these processes is likely to be negli gible com- 
pared to D production during BBN (but see also I JedamzikI 
12003). Stellar nuclear reactions cause a net destruction 



' ^ = 100 k,n /-iMpc-1 (H'lbble constant) 



of D, so measurements of the D abundance in sites that 
have undergone star formation, such as our Galaxy, pro- 
vide a lower limit to the priordial D/H value. However, 
many QSO absorption clouds are thought to have under- 
gone very little star formation, and we expect D/H in 
these clouds to be close to the primordial value. Unfor- 
tunately, very few absorption systems have the high col- 
umn density (> lO'^^ absorbers per cm^), simple velocity 
structure and narrow linewidth that are required for the 
Dl line to be separated from the nea rby Hi absorption 
JWebb. Carswell. Irwin fc Penstonlll99ll) . 

Most existing D/H measurements have been made in 
absorption systems with a redshift large enough that the 
Lyman limit falls in the visible {z > 2.2). For lower redshift 
systems the Lyman limit falls in the UV and must be ob- 
served from space. However, it is desirable to look for low 
redshift systems as there is a much larger sample of bright 
QSOs at lower redshifts. In addition, the number density of 
Lya forest absorbers is much lower at low redshifts com- 
pared to high redshifts. This means it is easier to accurately 
measure the continuum level around the relevant absorption 
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lines, and the chance of a randomly placed H I line contam- 
inating the D I absorption is smaller. 

The absorption system at z = 0.701 towards 
PG1718+4807 was first identified as an excellent can- 
didate for a D/H measurement from an International 
Ultraviolet Explorer (lUE) spectrum covering its Lyman 
limit (LL). Its LL is 'grey', meaning the flux does not drop 
to zero bluewards of the LL. The drop in flux at the grey 
LL strongly constrains the column density of the absorber. 
The column density was high enough for the D I Lya line to 
be detected, and the sharp break at the LL suggested the 
velocity strucure of the absorption complex was very simple. 



important for future analyses of D/H in QSO spectra. Fi- 
nally, we assess the likelihood that this absorption system 
shows D I and give an estimate of D/H for this system. 



2 DATA AND ANALYSIS 

The QSO PG1718-I-4807 has been observed with the Faint 
Object Spectrograph (FOS), Goddard High Resolution 
Spectrograph (GHRS) and STIS on th e HST, and the lUE 
iLanzetta. Turnshek fc Sandovallll993l) ^. We are interested 



ilLfinding the parameters of the H I and D I absorption. The 

IWebb CarsweU. Lanzetta. Ferlet. Lemoine. Vidal-Madiar fc Bowe4»i-i<^;^^° of the FOS and lUE spectra are much lower 



Jl997l) published the first analysis of D/H in the 
z = 0.701 system, using a GHRS spectrum covering 
the Lya and Sim lines a nd th e lUE spectrum cov- 
ering the LL. IWebb et all Jl997l. .1998) fitted the Hi 
absorption with a single cloud, assuming the redshift 
of the Sim line was the same as the redshift of the 
Hi absorption. They found D/H of 1.8 - 3.1 x 10"*, 
about ten times higher than D/H measured in other, 
higher redshift absorption systems (~ 2 — 4 x 10~^, see 
O'Meara. Tvtler. Kirkm an^ ^uzuki^ Pro chaska. Lubin &: Wolfd 
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tTytler, Burks, Lu. Fan. Wolfe fc S avage (1993.) 
quently published a further analysis of D/H in the same 
absorber, adding Keck/HIRES spectra of Mgll lines associ- 
ated with the z = 0.701 absorber. Using the Mgll lines they 
considered several models with different constraints on the 
Hi component's redshift. They found a larger range of D/H: 
8- 57 X 10"^ 



iKirkman . Tvtler. O'Meara. Buries. Lubin. Suzuki. Carswell 
EoOl ) (hereafter KTOB) published new HST STIS (Space 
Telescope Imaging Spectrograph) spectra covering the LL, 
Lya and Sim lines of the absorber. The spectrum of the 
LL has a resolution of 10 kms~^, and all the Lyman series 
lines from Lyman-e to Lyman- 16 were resolved (see Fig.0. 
This spectrum confirmed the presence of a single strong 
Hi component an d allowed its r e dshift to be measured 
precisely. However, iKirkman et alj i200lf) found additional 
weak absorption to the red of the main Hi component. This 
absorption was too weak to be seen in the higher order 
Lyman series, but was seen in the new STIS spectrum of the 
Lya line. They also claimed that the b parameter [b = \/2o", 
where a is the velocity dispersion) of the putative D I was 
larger than that expected for a gas with a temperature 
and turbulence calculated from the Hi and metal lines' 
b parameters. Based on these t wo points they cla imed 
that the absorptio n fitted as Dl in lWebb et alJ |^23) a^nd 
iTvtler et al.l lll999fl is contaminated by weak H I absorption. 
Due to this contamination they concluded the system does 
not provide evidence for a high D/H value. 

We present a further analysis of D/H in the z = 0.701 
absorption system, in which we look at the effect of system- 
atic errors, not addressed in previous analyses, that affect 
the z = 0.701 absorption system parameters. In particular 
we analyse how the the putative D I absorption's param- 
eters are affected by wavelength calibration errors in the 
HST spectra and present revised parameters for this partic- 
ular absorption spectrum. These calibration effects will be 



ra are 

(FWHM ~ 150 and ~ 900 kms"^ respectively) than those 
of the GHRS and STIS spectra. The FOS and lUE spectra 
do not provide any useful constraints on the parameters of 
D I or any H I sub-components which may be present in the 
z — 0.701 absorption complex. The parameters of the main 
H I component are very well determined by the STIS spectra 
alone (see section lTSll . Therefore we use only the GHRS and 
STI S spectra in our ana lysis. 

iTvtler et all (Il999l) detected Mgii with Keck/HIRES. 
We use their best fitting parameters for Mgll, with a mod- 
ified error for the Mgll b parameter (see section 17!^ . 



2.1 GHRS data 

The GHRS spectra were taken using the G140M grating 
with a 0.2 x 0.2 aperture giving a line spread function (LSF) 
with a FWHM of 14 kms~^. The observations consisted of 
56 exposures, which were shifted and added using the IRAF 
routines written for this purpose by the Space Telescope Sci- 

¥if& jiP^^ ^^VSrti&l^ ^^) ' P'^ff^^^^ ^^'^ specalign. We used the 
"pipeline reducf3d spectra, processed with CALHRS version 
1.3.13. The wavelength range covered is 2049 - 2098 A. The 



S/N per pixel near the Lya and Sim lines in the combined 
spectrum is ~10. 



2.2 STIS data 

The STIS observations were taken with the E230M echelle 
grating and the G140M longslit grating with the 0.2x0.2 
aperture. The FWHM of each of the LSFs of the STIS obser- 
vations are 2 pixels for the E230M echelle grating (5 kms~^ 
per pixel) and 1.5 pixels for the G140M grating (0.06 A or 
~11 kms~^ per pixel). Four exposures were taken with the 
G140M grating, and five with the E230M grating. In each 
case we used the pipeline reduced spectra (CALSTIS ver- 
sion 2.8). For each set of observations the exposures were 
rebinned to a common wavelength scale and combined us- 
ing a variance weighted average. The E230M spectra cover 
the wavelength region 1842 - 2673 A. The final S/N per 
pixel near the Lya and Sim lines in the averaged E230M 
spectrum is ~5. The G140M spectra cover the Lyman-e 
transition down to the Lyman limit, a wavelength range of 
1540 - 1594 A. The S/N per pixel in the G140M averaged 
spectrum is ~9. 



2 The GHRS, STIS, FOS and lUE spectra are all available f rom 
the HST multimission data archive: |http: / / archive. stsci.edu| 
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Figure 1. The STIS spectrum of the z = 0.701 Lyman hmit, a weighted sum of four STIS exposures with 
the G140M grating. There appears to be Galactic Civ absorption and emission at ~1548 and ~1550.5 A. 
We are unsure what is causing the absorption at ~1554.5, 1565 — 1566.5 and -^1580. 5 A. Neither these 
regions, nor the regions affected by Galctic Civ absorption were included in the fit. The normalised residuals 
(which we define = ^-^(data point — fitted value)/error) and the la error levels for the residuals are shown 
centred on y = —0.3. The thin curve shows the best fitting solution when the Lyo lines and the LL are fit 
simultaneously. This fit provides the H I parameters in Table 



2.3 Line spread functions 

The LSF of the GHRS observations can be approximated 
by a Gaussian, but the STIS LSFs are significantly differ- 
ent from a Gaussian. We used the instrumental profiles for 
the E230M and G140M gratings given in the STIS instru- 
ment handbook^. The actual STIS LSFs may be slightly 
asymmetric. The nature of this asymmetry depends on the 
orientation of the slit on the sky (Sahu in KTOB). We av- 
eraged the asymmetric sides of the LSF to give a symmetric 
LSF for both the G140M and E230M gratings. We note that 
even if some small asymmetry is present, it will not have a 
significant effect on the fitted absorption line parameters. 
Once the LSF per pixel for the original spectrum was deter- 
mined it was rebinned to the same wavelength scale as our 
final averaged spectrum. For the E230M grating, the pixels 
are a constant size in velocity rather than wavelength. We 
calculated an E230M LSF at the central wavelength of each 
line fitted in the E230M spectrum. 

2.4 Continuum placement 

The placement of the continuum in the STIS G140M spec- 
trum has a potentially significant effect on the fitted h pa- 

^ http: / /www. stsci.edu/hst/stis/performance/spectral_resolution/ 



rameter and column density of the main Hi component. To 
test the magnitude of this effect, we fitted the continuum of 
this region using several methods. Firstly we fitted a con- 
stant flux to regions apparently free of absorption. Secondly 
we fitted a power law fx oc A~". Finally, the lUE spectrum 
covers a larger wavelength range redwards of the Lyman 
limit than the STIS G140M spectrum, so it may give a bet- 
ter estimate of the continuum than just fitting the STIS 
Lyman limit. We fitted the lUE Lyman limit with a 3rd 
order Chebyshev polynomial and scaled this continuum to 
the level of the STIS Lyman limit, based on regions appar- 
ently free from absorption. Each of these methods gave very 
similar column densities and h parameters for the main H I 
component. The small error due to continuum fitting was 
added in quadrature to the statistical error. This combined 
error is the first error given in Tabled 



The error in continuum placement around the Lya line. 
Si III line and 4 doublet (Fig. |21 El and |1J was determined 
by fitting both straight lines and 3rd order Chebyshev poly- 
nomials to regions apparently free from absorption close to 
each line. The differences between these two methods of con- 
tinuum fitting have a negligible effect on the absorption line 
parameters. 
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Figure 2. The GHRS (above) and STIS (below) spectra covering the z = 0.701 Lyo line. The normalised 
residuals and their Itr error levels are shown centred on y = —0.4. The thin curve shows the best fitting 
solution when the Lya lines and the Lyman limit are fit simultaneously. This fit provides the parameters in 
Table HI The contributions from the main Hi component, red and blue components are shown by dashed 
curves. Note that the 'shape' of the STIS E230M Lyo line appears to be different to that of the GHRS line. 
In particular, the shape of the red wing and the absorption at ~ 2067.4 A are different. We are not certain 
what causes these differences, but they may be explained by correlations in the noise. 



2.5 Wavelength calibration and shifts between the 
HST spectra 

Misalignment of the wavelength scale between the GHRS, 
STIS G140M and STIS E230M spectra can potentially affect 
all the fitted parameters for H I and the putative D I, and the 
redshift for Sim and 4. 

The GHRS wavelength scale appears to be offset from 
the STIS wavelength scale by ~ 0.07 A (~ 1.4 GHRS pix- 
els). Despite helpful discussion with Glaus Leitherer from 
STScI we were unable to discover the reason for the discrep- 
ancy. KTOB shifted the GHRS spectrum to coincide with 
the STIS spectrum, finding the amount by which to shift 
by cross correlation over the common wavelengths covered. 
They then checked this shift by comparing the positions of 
profile fits to sharp features in each spectrum. However, they 
do not estimate the error in their shift. 

We calculate the shift by choosing the three sharpest 
absorption features at relatively high S/N present in each 
spectrum and calculating the cross correlation for each of 
them. The regions we used were: 2051.1 — 2053.7 A (the 
Sim line), 2056.5 - 2060.4 A (a Ly/3 line at z = 1.0065) 
and 2069.6 - 2072.2 A (a Lya line at ^ = 0.70326). The 
three shifts for these regions were 0.070 A, 0.072 A and 
0.058 A respectively. The three shifts are roughly consistent, 
suggesting that the wavelength offset is the same across the 



spectrum. We take the average of these three shifts and move 
the GHRS spectrum by this amount to align it with the STIS 
spectrum, assuming that the STIS wavelength calibration is 
more accurate than the the GHRS calibration. To quantify 
the error in this shift, we fit the sharpest feature in each 
spectrum, the Sim line, with VPFIT. The VPFIT errors in 
the position of each line were then added in quadrature to 
give an estimate of the total error in the shift. This error 
is ~ 0.007 A, or ~ 1 kms~^. This error does not a have a 
significant effect on the the parameters of the putative D I 
or main Hi absorption. 

For a 0.2x0.2 aperture and the MAMA detector used 
for the STIS observations, the absolute wavelength calibra- 
tion error can be as much as 0.5 - 1.0 pixels (2a. IProffitd 
ioo^). In our analysis we introduce a conservative error of 
0.3 pixels in the G140M and E230M wavelength calibrations. 
Thus we consider three different cases: (1) The wavelength 
cafibration of both the G140M and E230M spectra are both 
correct, (2) the E230M spectrum wavelength scale shifted 
0.3 pixels bluewards and the G140M wavelength scale 0.3 
pixels redwards, and (3) the E230M wavelength scale shifted 
0.3 pixels redwards and the G140M scale 0.3 pixels blue- 
wards. Introducing even these conservative relative shifts of 
0.6 pixels between the STIS E230M and G140M spectra does 
have a significant effect on the absorption line parameters, 
particularly those of the putative D I. This is because the 
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Figure 3. The GHRS (above) and STIS (below) spectra covering the z = 0.701 Sim hne. The normahsed 
residuals and their Icr error levels are shown centred on j; = 0. The thin curve shows the best fitting solution 
when both lines are fit simultaneously. This fit provides the parameters in Table [Tl 




T 



"T r- 



2632 



2634 



2636 



2638 



Wavelength (A) 



Figure 4. A tentative detection of the z = 0.701 4 doublet in the STIS E230M spectrum. The normalised 
residuals and their Icr error levels are shown centred on j/ = 0. The thin curve shows the best fitting solution 
when both lines of the doublet are fit simultaneously. This fit provides the parameters in Table [Tl 



redshift of the main H I component is determined by the 
higher order Lyman series lines in the G140M spectrum. In- 
troducing a shift between the G140M and E230M spectra 
means that the position of the main Hi absorption in the 
Lya line changes. To compensate for this change in position, 
the parameters of the putative D I and the red H I compo- 
nent are also changed. We describe the effect of these shifts 
on the parameters of the putative D I in section [ 



2.6 Contamination of the z = 0.701 system by 
systems at other redshifts 

We attempt to identify all lines in the E230M spectrum to 
determine if the z — 0.701 absorption lines are contaminated 
by H I or metal lines from systems at different redshifts. 
There are six absorption systems with log[A^(Hl)]> 14.5 
present in the E230M spectrum. They are at redshifts of 
1.0867, 1.0674, 1.0549, 1.0318, 0.7247 and 0.5272. The Ly/3 
line at 2: = 0.5272 falls at ~ 1566. 5A, and may contribute 
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Figure 5. A section of the STIS G140M spectrum showing the 
constraints on the column density of the z = 0.602 absorption 
system. The thin solid curve is the best fit obtained by including 
all Hi components from the z = 0.70108 system and a single 
component at z = 0.6023. The middle dashed curve shows the 
z = 0.6023 column density set at the 2(t upper limit. The lower 
dashed curve shows t he col umn density set to log[A'^{H l)]= 16.7, 
given inlrvtlcr et aT] il99Hh . The h parameter for the z = 0.6023 
line cannot be well determined by the data and was fixed at 15 
kms~^. Increasing the b parameter decreases the upper limit on 
the z = 0.6023 column density. 

to the absorption seen at this wavelength in the G140M 
spectrum (see Fig. Ql. We did not include the region with 
this absorption when we fitted the G140M spectrum. For 
each log[A^(Hl)]> 14.5 system we searched the wavelength 
ranges covered by the G140M and E230M spectra for strong 
associated metal lines that may contaminate the z = 0.701 
absorption lines, including Sill, Niv, Ol and Cll. Only two 
systems show any strong metal lines: the z = 1.0867 system 
has associated Sim absorption, and the z = 1.0318 system 
shows Nv absorption. Galactic Fen and Mnll absorption is 
also present. We found no other strong metal lines or higher 
order Lyman Hi transitions that may be blended with the 
z — 0.701 lines and affect our analysis. 

2.7 z = 0.602 absorption system 

iTvtler et al.l lll999ll suggest that there is evidence in the lUE 
spectra for another Lyman limit system with log[A''(Hl)] ~ 
16.7 at 2 = 0.602. It is important to determine whether 
such a system is present, since if it is, its Ly-7 line will fall 
at ~ 1558 A and may blend with the Ly-13 or Ly-14 lines 
of the z = 0.701 system. Such a Lyman limit system should 
have a Lya fine at ~ 1947 A. The STIS E230M spectra 
covers this region, although the signal to noise is poor (~ 3). 
A strong line does appear at 1948 A, but the low S/N and 
contamination from a higher redshift H I Ly/3 line mean the 
column density cannot be determined from this line alone. 



However, we can put an upper limit on the column den- 
sity of the suspected z = 0.602 system in the following way. 
We fitted the z — 0.701 Lyman limit and higher order lines, 
excluding the lines possibly blended with the z — 0.602 Ly-7 
line. Then we compared the predicted Ly-13 and Ly-14 lines 
with the absorption seen, assuming any residual absorption 
in these two lines is due to the z = 0.602 Ly-7 line. Com- 
bining the constraints from the low S/N Lya line and the 
method above, we find that the contaminating system has a 
redshift of 0.6023± 0.00003, placing it between the z = 0.701 
Ly-13 and Ly-14 lines with log[iV(Hl)]< 14.54 (99.7% confi- 
dence). The z = 0.602 system does not affect the measured 
parameters of the z = 0.701 system, as illustrated in Fig. 12. 71 
No metal lines were detected at z = 0.602. 



2.8 Voigt profile fitting 

Voigt profiles were fitted to the absorption lines and the best 
fitting parameters were determined by minimising the 
statistic. The program VPFIT'* was used for all line profile 
fitting. The absorption lines that were detected and fitted 
in the z = 0.701 system in the STIS spectra are the Lya 
line, the Lyman series from Ly-e to the Lyman limit, Sim 
(1206.5 A) and the 4 doublet (1548.2 and 1550.8 A). The 
GHRS spectrum covers the Lya and Sim lines. We fitted 
all parameters for case (1) described in section [2.51 for no 
shift between the STIS spectra. The 0.3 pixel error in the 
uncertainty in the STIS wavelength calibrations is listed af- 
ter the statistical error in Table 1 for the lines it applies to. 
We look at the effect of the shifts between the STIS spectra 
wavelength calibrations on the putative D I's parameters in 
section 12.91 

The parameters of the main Hi component are tightly 
constrained by the unsaturated H I higher order lines and 
the drop in flux at the grey Lyman limit. The main uncer- 
tainties in this case are the position of the continuum and 
the accuracy of the STIS wavelength calibration. Table Q 
shows the best fltted parameters to the main Hi compo- 
nent. The error in redshift due to the error in STIS zero 
point wavelength calibration is given after the statistical er- 
ror. The first error for each parameter given in the table is 
the quadrature addition of the statistical fitting error and 
the continuum placement error, where applicable. 

The Sim line parameters were found by fitting the 
GHRS and STIS spectra simultaneously. Again we include 
separately the statistical error in redshift and the error 
when the wavelength calibration uncertainty of the STIS 
E230M spectrum is taken into account. The Sim line is off- 
set bluewards from the main H I system by 0.5±1 (statistical) 
kms~^. 

The 4 fine parameters were found by fltting the 4 dou- 
blet in the STIS E230M spectra. The statistical error and 
error due to wavelength calibration are shown separately. 
The 4 lines are offset bluewards from the main Hi compo- 
nent by 4 ± 3 (statistical) kms~^. The best fitting h pa- 
rameter for 4 is large compared to 6(Mgll) and ti(Silll). 
This difference is also seen in many damped Lya absorbers 
JWolfe fc Prochaskallioool) and is probably due to 4 having 
a much higher ionization energy than Sim and Mgll, and 

^ Carswell et al., http//www. ast.cam.ac.uk/~rfc/vpfit. html 
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Table 1. The absorption line parameters and their errors given by VPFIT. For the two Hi and putative Dl 
lines these are the best fitting parameters when the STIS Lyman limit abd the GHRS and STIS Ly« lines 
are fit simultaneously. For Si III these are the best fitting parameters when the GHRS and STIS Si III linos 
are fitted simultaneously. The first errors given are the quadrature addition of the statistical and continuum 
fitting errors. The second errors for the Sim line, 4 doublet and main Hi component are the errors due 
to uncertainty in the wavelength calibration of the STIS spectra (0.3 pixels). Note that the errors for the 
putative D I and H I red sub-component are not meaningful (see section 12.91 . We list the Mg II parameters 
published by KTOB based on observations taken with Keck/HIRES. The error for fe(Mgll) is our estimate 
of the error based on our simulated spectra (see section f2.8l . The error for 6(Mgll) in KTOB is 1.1 kms~^. 



Ion 


log(Ar) (cm-2) 


z 


b ( kms 


Hi 


17.213 ± 0.007 


0.701074 ± 0.000003 ± 0.000015 


22.1 ± 0.4 


Hi " 


13.5 


0.70142 


19 


Putative D I ^ 


13.88 


0.701108 


22 


Si III 


12.84 ± 0.05 


0.701071 ± 0.000007 ± 0.000009 


14 ± 2 


4 


13.18 ± 0.08 


0.701051 ± 0.000019 ± 0.000009 


20 ± 5 


Mgii 


11.50 ± 0.05 


0.701100 ± 0.000005 


12.0 ± 1.4 



VPFIT does not return sensible errors for this H I component. The reasons for this are described in section 
^ See section [2.9l Fig. ISllTl andlSlfor the putative Dl parameter errors. 



so being associated with a different velocity space in the 
absorption cloud. 

We use the parameters of the Mg ll lines given in KTOB. 
To verify the errors given by KTOB, we generated a syn- 
thetic Mgll spectrum from the fitted b, N and z values 
they provided, at the same S/N (~70) and resolution (8 
kms~^ FWHM) as the Mgll spectrum they published. Our 
a-{N) and (j{z) were very similar to KTOB's, but our a{b) 
was larger, 1.4 kms~^ instead of 1.1 kms~^. We use our 
larger a{b) estimate in our analysis. The Mgll lines are off- 
set redwards from the redshift of the main Hi system by 
5± 1 (statistical) kms~^. The systematic error for the wave- 
length calibration for Keck/HIRES spectra is < 0.5 kms~^. 
This 5 kms~^ shift may be due in part to an error in the 
STIS wavelength calibration (the 0.3 pixel error in wave- 
length calibration we consider corresponds to a velocity er- 
ror of ~3 kms~^). Even though we are not sure of the cause 
of this shift, it does not affect our analysis, since we use only 
the Mgll 6 parameter (section l2.10ll . 

We attempted to fit the LL and both Lya lines with 
a single Hi component. The minimum reduced = 
Xmin/(iio. of degrees of freedom) for this fit was 1.7. This 
fit predicted too much absorption in the higher order lines, 
and could not account for extra absorption in the red and 
blue wings of the Lya line. We needed to add two extra ab- 
sorption lines, one H l sub-component to the red of the main 
component, and a putative D I line to the blue, to find an 
acceptable reduced (1.04). Table Ogives the parameters 
of the red and blue H I components. 



2.9 The errors on the parameters for the putative 
Dl 

The la errors on the parameters in Table 1 are those calcu- 
lated by VPFIT. These errors are based on the assumptions 
that (1): the model we have chosen (a single main Hi com- 
ponent with a red and blue sub-component on either side) is 
correct and (2): the fitted parameters are independent and 



the shape of the x^ parameter space around the minumum 
X^ is parabolic. However, for the putative D I the column 
density, b parameter and redshift are correlated. In addi- 
tion, the x^ parameter space is asymmetric (see Fig. |S1 
and (HI . This is due to the small number of data points that 
are available to constrain the parameters of the putative D I. 
It may also be due to the difference in 'shape' between the 
STIS and GHRS Lya lines (see Figurej^J. For these reasons, 
the la estimates given by VPFIT for the putative D I pa- 
rameters may not be accurate. This is also likely to be the 
case for the red Hi sub-component. For the other lines (the 
main Hi component, Sim, 4, and Mgll), the la errors given 
by VPFIT will be sufficiently accurate. 



We find a robust parameter error estimate for the D I 
parameters by using Ax^ probability contour maps. These 
are generated by fixing two parameters (for instance, the b 
parameter and redshift of the putative D l) and then vary- 
ing the remaining parameters (the column density of the 
putative Dl and all the parameters of the main Hi compo- 
nent and red Hi component) to minimise x^- This creates 
a grid of x^ values. The distribution of Ax^ = ~ Xmin 
is the same as that of a x^ distribution with a number of 
degrees of freedom equal to the number of fixed parameters 
(in this case, two). Here Xmin is the smallest value of x^ for a 
particular grid. Provided that errors on the data points are 
normally distributed, the probability contours generated in 
this way will be correc t, irrespect ive of correlations between 
the fitted parameters iCashlll97(J) . We also need to consider 
the effect the shifts between the STIS E230M spectra have 
on the D I parameters. To do this we again consider the 
three shifts described in the previous section. Fig.|S]|7|and|5] 
show the probability contour maps for the three parameters 
of the putative D i for the three different relative shifts be- 
tween the STIS Lyman limit and Lya spectra, representing 
the uncertainty in the wavelength calibration between the 
two spectra. 
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Figure 6. logio(Af[H l]) versus velocity for the putative Dl. For the left graph the Lyman limit is shifted 
0.3 pixels redwards and the Lyo lines are shifted 0.3 pixels bluewards. For the right graph the Lyman limit 
is shifted 0.3 pixels bluewards and the Lya lines are shifted 0.3 pixels redwards. For the centre graph the 
Lyman limit and LyQ lines are fitted with no shift. The velocity zero point corresponds to 81.6 kms~^ 
blueward of the best fitting position of the main H I component, where we expect D I to be. The contours 
represent, from innermost to outermost, the 68.4%, 95.4%, 99.73% and 99.99% confidence levels. The vertical 
hashed region represents the 68.4% confidence region on the positon of the main Hi component. 
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Figure 7. b parameter versus velocity for the putative D I. For the left graph the Lyman limit is shifted 0.1 
pixels redwards and the Lya lines are shifted 0.3 pixels bluewards. For the right graph the Lyman limit is 
shifted 0.3 pixels bluewards and the Lya lines are shifted 0.3 pixels redwards. For the centre graph the Lyman 
limit and Lya lines are fitted with no shift. The velocity zero point corresponds to 81.6 kms~^ blueward of 
the best fitting position of the main H I component, where we expect D I to be. The contours represent, from 
innermost to outermost, the 68.4%, 95.4%, 99.73% and 99.99% confidence levels. The vertical hashed region 
represents the 68.4% confidence region on the positon of the main Hi component. The horizontal hashed 
region represents the 68.4% confidence region for the D I b parameter predicted by the H I, Si III and Mg II 
(see Fig.inj. 



2.10 The putative D I's b parameter 

We can predict what the 6 parameter of a D I Une associated 
with the main Hi absorption from the b parameters of the 
Sim, Mgll and main Hi component. This is done by mea- 
suring the thermal line broadening, bthcrm, and turbulent 
broadening, 6turb, by using the b parameters from all avail- 
able ions. If we assume the absorbing cloud has a thermal 
Maxwell-Boltzmann distribution and any turbulence can be 
described by a Gaussian velocity distribution, then the 6 
parameter for a particular ion will be given by 

bion ~ ^thcrm + ^turb- (1) 

Here &therm = where T is the temperature of the 
gas cloud, m is the mass of the absorbing ion and k is Boltz- 
mann's constant, fetm-b represents the Gaussian broadening 



due to small scale turbulence, and is the same for all ionic 
species. &?hcrm is proportional to the inverse of the ion mass. 
We plot b^ against inverse ion mass in Fig. |^ Subject to 
the assumptions above, all the ions in the same cloud veloc- 
ity space should lie on a straight line whose intercept gives 
&turb and slope gives the cloud temperature. Since Mgll and 
Hi have very similar ionization potentials (IP) (13.6 eV and 
15.0 eV respectively), they should trace out the same gas. 
The IP of Sim is somewhat higher (33.5 eV), but there is 
evidence that 'intermediate' IP ions, such as Aim, have a 
similar velocity structure to lower IP ions, such as Mgll 
J Wolfe fc Prochaskallioool) . We assume that the Mgll and 
Si III hnes are in the same gas as the H I and fit least squares 
line of best fit to these three points. This is shown as the 
dashed line in Fig. |n| If we use KTOB's estimate for Mgll 
a{b) . The D I point is plotted with its la error bars as given 
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log(N[DI]) log(N[DI]) log(N[DI]) 

Figure 8. b parameter versus logio (-'^P l]) for the putative D l. For the left graph the Lyman hmit is shifted 
0.3 pixels redwards and the Lya hnes are shifted 0.3 pixels bluewards. For the right graph the Lyman hmit 
is shifted 0.3 pixels bluewards and the Lya lines are shifted 0.3 pixels redwards. For the centre graph the 
Lyman limit and Lya lines are fitted with no shift. The contours represent, from innermost to outermost, 
the 68.4%, 95.4%, 99.73% and 99.99% confidence levels. The horizontal hashed region represents the 68.4% 
confidence region for the Dl b parameter predicted by the b parameters of Hi, Sim and Mgll (see Figl9l. 
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Figure 9. The b parameter squared ( kms~^)'^ versus the inverse 
ion mass (amu~^). The dashed line is the least squares line of best 
fit to the Mgll, Sim, and Hi points. 4 and the putative Dl are 
shown for comparison. The error bars in each case are calculated 
from the Icr errors given by VPFIT. Note that the D I error bars 
are not particularly meaningful (see section 12.91 . The contour 
plots shown in Fig ISl 171 and ISl are required to determine whether 
f)(Dl) is consistent with f)(Hl), fe(Mgll) and f)(Siin). 



2.11 Is the blue Hi component deuterium? 

If the blue component is D I, it should fall at a redshift cor- 
responding to a velocity difference of 81.6 kms~^ bluewards 
of the main Hi component. In addition, it should have a b 
parameter consistent with that predicted in the last section. 

We use the contour plots described in section TIM com- 
bined with the b parameter test in the last section, to see if 
the putative Dl component's parameters are consistent with 
those of D I. Fig. |n] Q a^nd |5] show how A^^ changes for a 
range of b parameter, redshift and column density values. 
Note the degeneracy between the parameters, which has a 
consequence that the errors obtained using VPFIT are not 
likely to be accurate for the Dl. The hashed bands on each 
of the plots show the 68% confidence region that the param- 
eters of the blue component are expected to fall in if it is D I. 
Fig.|7|best illustrates the situation — if the blue component 
is D I, we would expect its redshift to be the same as that of 
the main Hi component (0 kms~^ velocity difference) and 
the b parameter to be consistent with the b parameters of 
the Hi and metal lines (~ 17 kms"'^). 

For the left case of Fig. |7| corresponding to a fixed rel- 
ative shifts between the G140M and E230M spectra, we can 
see that the b parameter is completely consistent with that 
of D I, but the redshift is inconsistent. For the right case, 
the redshift is consistent, but the b parameter is not. Where 
there is no relative wavelength shift between the spectra, 
the parameters of the blue component are inconsistent at 
a ~ 80% confidence level. While these results do not sug- 
gest that the blue component is D I, they do not conclusively 
show that is is not D I. 



by VPFIT. We find a temperature of 2.05 ±0.12 x lO'^K and 
a turbulent broadening of 12.1±1.3 kms~^. Although we use 
our estimate of Mgll a{b) rather than KTOB's, if we do use 
KTOB's Mgii (t(6), the derived temperature and turbulent 
broadening do not change significantly. 



If we take the case of no shift between the G140M and 
E230M spectra and assume the blue absorption is entirely 
due to D I, fixing the b parameter and redshift at those we 
expect for Dl, we find the ratio D/H = (3.0 - 4.6) x 10"" 
(la limit). 
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2.12 Comparison with previous analyses of this 
system 

We have analysed this absorption sys t em using the same 
spectra presented in lWebb etall ^ll997^ . lTvtler et alj (Il999l) 
and KTOB, considering important systematic effects not ad- 
dressed in these previous analyses. 

There are several small differences in our fitted param- 
eters for the main Hi component, putative Dl component 
and Sim line compared to KTOB's fitted parameters. Our 
line spread function is slightly different from KTOB's and 
we use a larger error for &(Mgll) (see section IT^ . We have 
explored the effect that relative shifts between the wave- 
length scales of the G140M and E230M STIS spectra have 
on the parameters of the putative D I line. However, none of 
these effects substantially change the conclusions of KTOB. 
As can be seen from Fig. |3 if a large relative shift between 
the G140M and E230M wavelength scales is present, the 
probability that the putative D I really is D I is decreased. 

The most important reason for our different result from 
KTOB is a difference in generating our probability contours 
in Fig.El|3and|Hl When generating our x contours, we var- 
ied all the parameters of all the components in the z = 0.7 
complex. To generate the plot in KTOB's Fig. 7, when cal- 
culating the minimum only the column density of the 
putative D l component was varied. The remaining putative 
D I parameters and all the parameters of the main H I com- 
ponent and the red H I sub-component were fixed by KTOB. 
This means their contours do not accurately represent the 
actual range of probabilities for the putative D I b parame- 
ter and redshift. The true confidence ranges are somewhat 
larger. This is the main reason we find that the putative D I 
parameters are inconsistent with those expected for D I at 
the ~ 80% rather than > 98% level. 

Our D/H upper limit calculated assuming the blue ab- 
sorption is entirely due to D I is, not surprisingly, consistent 
with the lTvtler et al. ( 1999[) upper limit of 5.7 x 10"'*. 



3 DISCUSSION 

We have analysed the z = 0.70f absorption complex to- 
wards PG 1718-H4801 using HST STIS and GHRS spectra, 
considering systematic errors that were neglected in previ- 
ous analyses. The most important of these systematic errors 
is the absolute wavelength calibration of the STIS Gf40M 
and E230M spectra. Introducing a relative shift, comparable 
to the wavelength error, between these two sets of spectra 
significantly alters the parameters of the putative D I. 

We find that the parameters of the absorption line pre- 
viously identified as D I are marginally inconsistent with 
those expected. In particular, the 6-parameter and line posi- 
tion of this feature are at best consistent with the expected 
values at the 20% level. 

We have shown that the absorption system at z = 0.602, 
previously claimed as having log[A'^(Hl)] ~ 16.7 in fact has 
log[Ar(Hl)] < 14.5. 

If we assume that the blue component is D I, we find 
D/H = (3.0 — 4.2) X 10"*. If significant contamination is 
present, then of course D/H may be significantly lower than 
3.0 X 10"*. 

The D/H range given here corresponds to a Icr range 



for Q,bh? of 0.003 — 0.006. We can compare this range with 
Q,bh^ estimates from other sources. Two recent estimates of 
the primo rdial Li^/H abundance arc (0.91 — 1.91) x 10"^" 
iRva n. Beers. Olive. Fields fc Norr is 2000) and (1.92 - 
2.49) X 10" ^° llBonifacidl2002l) . Combining these two ranges, 
we find U'' {Vtbh^) = 0.0044 - 0.020. Two samples recently 
used to measur e the He"^ abundance, Y p, give 0.238 ± 0.003 
(The sample in 'izotov & Thuan"'l998', using the ionization 
correction factor iii,GnMnwal^_^teigmM^^^/ieKas 200*3) 
and 0.2405±0.0017 jPeimbert. Peimbert fc Ruizl 2000'). Tak- 
ing the highest and lowest limits of these two ranges gives 
He*(nb/i^) = 0.0065 - 0.0145. Sievers et al. (2002) give 
CMB{nth'^) = 0.023 ± 0.003 based on aU currently avail- 
able CMB data. Taking an average of the five D/H mea- 
surements made in other QSO absorption clouds, weighted 
by their inverse variances, gives D{Qbh^) = 0.021 ± 0.001. 

The Qbh^ range for the z = 0.701 QSO absorber is 
consistent with the Li'^ range and only marginally incon- 
sistent with the He^ range. It is significantly inconsistent 
with both the CMB and low QSO D/H Qbh^ ranges, how- 
ever. The mctallicity of this absorber is low, [Si/H] — —2.4 
(Tvtlcr ct al. 1999), and the QSO absorption systems used 
to measure D/H have similarly low metallicities. At these 
metallicities no significant astration is thought to have oc- 
curred and the measured D/H values should be primordial. 
For a standard homogenous BBN, we expect the primor- 
dial D/H to be independent of direction and redshift. Thus 
the inconsistency of the z = 0.701 D/H value with lower 
D/H values in other QSOs and the CMB Qbh^ value sup- 
ports the argument that the putative Dl in this system is 
contaminated with Hi. 

However there are several caveats. Only five D/H mea- 
surements in QSO absorbers have been made. Even among 
these 'low' D/H values there is considerable scatter, which 
is not expected based on the standard BBN, no astra- 
tion picture. This is may be due to unaccounted for sys- 
tematic errors, but other explantions have been s uggested. 
iFields. Olive. Silk. Casse fc Vangioni-FlamI fcOOll) suggest 
the scatter may be evidence for an early population of stars 
that ha ve caused a s tratio n even in these low metallicity ab- 
sorbers. Ijedamzi^ 1120021) suggest it is possible that there 
may be rare astrophysical sites where D/H has been en- 
hanced above the primordial value. It would be imprudent to 
reject entirely the possibility that this system shows a high 
D/H until further D/H measurements in QSO absorbers are 
available, particularly if there is a trend of D /H with metal- 
licity. 

We acknowledge helpful comments from Michael Mur- 
phy and Mike Irwin, and correspondence regarding STIS 
wavelength calibration with Glaus Leitherer at the STScI. 
We also thank Gary Steigman for his helpful correspon- 
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the paper. 



REFERENCES 

Adams T. F., 1976, A&A, 50, 461 
Bonifacio P., 2002, A&A, 395, 515 
Cash W., 1976, A&A, 52, 307 

Epstein R. I., Lattimer J. M., Schramm D. N., 1974, Nat, 
263, 198 



The Deuterium Abundance in the t. — 0.7 absorber towards QSO PG1718+4807 



Fields B. D., Olive K. A., Silk J., Cassc M., Vangioni-Flam 

E., 2001, ApJ, 563, 653 
Gruenwald R., Steigman G., Viegas S. M., 2002, ApJ, 567, 

931 

Izotov Y. I., Thuan T. X., 1998, ApJ, 500, 188 

Jedamzik K., 2002, Planetary and Space Science, 50, 1239 

Kirkman D., Tytler D., O'Meara J. M., Buries S., Lubin 

D. , Suzuki N., Carswell R. F., Turner M. S., Wampler 

E. J., 2001, ApJ, 559, 23 

Kirkman D., Tytlcr D., Suzuki N., O'Meara J., Lubin D., 

2003, astro-pli/0302006 
Lanzetta K. M., Turnshek D. A., Sandoval J., 1993, ApJS, 

84, 109 

Lcvshakov S. A., Dcssaugcs-Zavadsky M., D'Odorico S., 

Molaro P., 2002, ApJ, 565, 696 
O'Meara J. M., Tytler D., Kirkman D., Suzuki N., 

Prochaska J. X., Lubin D., Wolfe A. M., 2001, ApJ, 552, 

718 

Peimbert M., Poimbcrt A., Ruiz M. T., 2000, ApJ, 541, 688 

Pettini M., Bowen D. V., 2001, ApJ, 560, 41 

Proffitt C, 2002, STIS Instrument Handbook, Version 6.0. 
STScI, Baltimore STScI 

Ryan S. G., Boers T. C., Olive K. A., Fields B. D., Norris 
J. E., 2000, ApJ, 530, L57 

Tytler D., Buries S., Lu L., Fan X., Wolfe A., Savage B. D., 
1999, AJ, 117, 63 

Webb J. K., Carswell R. F., Irwin M. J., Penston M. V., 
1991, MNRAS, 250, 657 

Webb J. K., Carswell R. F., Lanzetta K. M., Ferlet R., 
Lemoine M., Vidal-Madjar A., 1998, in Structure et 
Evolution du Milieu Inter-Galactique Revele par Raies 
D'Absorption dans le Spectre des Quasars, 13th CoUoque 
d'Astrophysiquc dc I'lnstitut d'Astropliysiquc dc Paris A 
High Deuterium Abundance at z=0.7; Evidence for Cos- 
mic Inhomogeneity?. pp 345 — h 

Wcbb J. K., Carswell R. F., Lanzetta K. M., Fcrlct R., 
Lcmoinc M., Vidal-Madjar A., Bowcn D. V., 1997, Nat, 
388, 250 

Wolfe A. M., Prochaska J. X., 2000, ApJ, 545, 603 



